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This work explores the scaling effects for FPSE (free piston Stirling engines), which are known for their 
simple architecture and potentially high thermodynamic performances. Scaling laws are given and their 
potential for miniaturization is highlighted. A simple model which allows the design of the geometrical 
parameters of the heat exchangers, the regenerator and the masses of the pistons is proposed. It is based 
on the definition of six characteristic dimensionless groups. They are derived from the physics underlying 
the behavior of the free piston Stirling machine and their relevancy is backed up by comparisons be¬ 
tween documented Stirling engines from the literature. Keeping constant values for each group 
throughout the scaling range theoretically ensures constant performance. The main losses of Stirling 
engine (heat conduction loss, reheat loss in the regenerator, pressure drop and gas-spring hysteresis) can 
be expressed as a function of the geometrical and operating parameters. Additionally, the consequences 
of leakage due to the manufacturing precision of pistons architectures are underlined. From the proposed 
scaling laws, potential power and efficiency of Stirling cycle engines at a millimeter scale can be antic¬ 
ipated. It appears that the power density increases with miniaturization. It is also shown that the dy¬ 
namic masses related to the engine size are increased when scaling down and that the gap leakage 
presents the highest detrimental effects on the efficiency. These results call for dedicated architectures 
for micro-engines. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

The miniaturization of thermo-mechanical devices has become 
an increasing trend to match the portable power demand. With 
moderate efficiencies, micro-engines using hydrocarbon fuels 
combustion as a heat source make the realization of high density 
powering system possible [1—3]. They aim at delivering electrical 
power ranging from dozens of milliwatts to a few watts. These 
devices could be manufactured by MEMS (microelectromechanical 
systems) process approaches and then be ultimately integrated to 
global units of sensors, actuators and microprocessor as “power 
plants on a chip”. 

The Stirling cycle is one of the most interesting choices because 
of external heat source and no valves. It offers different possible 
applications: engine, cooler or heat pump. The FPSE (free piston 
Stirling engine) architecture is favorable for miniaturization 
because of its simple mechanical arrangement and the lack of 
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complex kinematic linkage. Fig. la presents a generic scheme of a 
FPSE which underlines the basic components: 

• Two heat exchangers (heater, cooler). They allow the gas to be 
kept at constant high 7h and low T c temperatures. 

• The regenerator. This component is usually made of a porous 
matrix. It periodically stores and releases heat from and to the 
gas, lowering the external heat needed to keep the hot gas 
temperature. 

• Two pistons (classically denoted piston and displacer). They 
ensure the proper dynamics of the FPSE. The displacer transfers 
the gas to the cold or the hot side of the engine (see arrows in 
Fig. 1) whereas the mechanical work is extracted from the 
piston motion. Note that the displacer also ensures a thermal 
isolation between the hot and cold sides of the engine. 

Based on a given Stirling engine, models have been proposed for 
optimization purposes [4-7]. The complexity is increased in the 
case of FPSE [8,9] because of the strong coupling between the dy¬ 
namic and thermodynamics equilibriums. The effect of a drastic 
miniaturization on the performances cannot be easily performed 
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Nomenclature 

Vi 

volume of component i (m 3 ) 



^sw 

displacement (m 3 ) 


displacer area (m 2 ) 

^swe 

expansion space swept volume (m 3 ) 


piston area (m 2 ) 

^swc 

compression space swept volume (m 3 ) 

-^rod 

rod area (m 2 ) 



Aff 

free flow area (m 2 ) 

Greek symbols 

c d 

displacer mechanical dissipation (N m -1 s) 

a 

swept volume phase angle (rad) 

C P 

piston mechanical dissipation (N m -1 s) 

A p 

Pressure loss (N m“ 2 ) 

d w 

wire diameter of the matrix regenerator (m) 

e 

scale (-) 

D d 

displacer diameter (m) 

Vr 

regenerator efficiency (-) 

f 

frequency (Hz) 

K 

swept volume ratio (-) 

/f 

friction factor (-) 

P 

dynamic viscosity (N m -2 s) 

h 

convection coefficient (W m -2 K -1 ) 

P 

material density (kg m 3 ) 

k 

thermal conductivity of gas (W m 1 K -1 ) 

CO 

pulsation (rad s -1 ) 

Li 

length of component exchanger (m) 

Ci 

dimensionless mass in component i (07 = m^M) 

f'ref 

reference length (L ref = U sw 1/3 ) (m) 


Dimensionless mass flow (cr f = rhj/(coM)) 

M 

total mass of gas (kg) 

1 

porosity of the regenerator (-) 

m d 

mass of a displacer (kg) 



m p 

mass of a piston (kg) 

Subscripts 

rrii 

mass of a gas within component i (kg) 

bp 

piston gas spring space 

n x 

number of channels (-) 

bd 

displacer gas spring space 

1 m 

mechanical power (W) 

c 

cooler 

P 

pressure (N m -2 ) 

d 

displacer 

R 

ideal mass gas constant (J kg -1 K -1 ) 

h 

heater 

r h i 

hydraulic radius of exchanger (m) 

P 

piston 

T c 

cold temperature (K) 

r 

regenerator 

T h 

hot temperature (K) 

ref 

reference 

T w 

wall temperature (K) 

rod 

Rod of displacer 

Uj 

mean gas velocity in component i (m s -1 ) 




using these models. Then, a dedicated approach is required to 
obtain design scaling rules and provide relevant information for 
novel architectures. 

Only a few works dealt with the specific issue of the size effect 
on the performances of miniature Stirling machines. Nakajima [10] 

a) 


performed a scale analysis on each of the elementary effects which 
govern the behavior of the Stirling engine but these effects were 
considered independently. 

Many works focused on the study of the regenerator because of 
the drastic repercussion of its effectiveness on the efficiency of the 

b) 


Gas spring for 
displacer 
(V bd ) 


Cooler 

(j'ho L c , n c ) 


Regenerator 
(r hn L n 72 ,) 


Heater 

(t hh T/j, 77/,) 



I 


Gas spring for piston 
(V bp ) 


Gap between 
the displacer and 
the casing 


Piston 

{tn p ) 


Displacer 

(m d ) 




B-B 

regenerator 


77 /, 



A-A 

heater 


Fig. 1. a) Simplified geometry for the study, b) main parameters. 
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Stirling machine. For cryocoolers, Peterson [11] developed an 
approach of scaling effects on the regenerator effectiveness and the 
axial thermal conduction. Yang and Chung [12] extended the study 
of the regenerator to take into account the void volumes, as well as 
the frictional and the viscous losses. Conclusions suggest that 
millimeter size machines are viable providing that low conductivity 
material is used. Nonetheless, because the considered pressure and 
the frequency were kept constant in these studies, restrictive de¬ 
ductions may have been obtained. The operating parameters 
(pressure and frequency) are indeed likely to be modified during 
scaling. 

Simple modeling approaches that can be performed without 
computational needs are called FOM (first-order methods) ac¬ 
cording to the Martini’s nomenclature [13]. They represent effective 
and useful models for preliminary design. The simplest FOM gives a 
relationship between the engine output power the mean pressure, 
the operating frequency, the displacement and a dimensionless 
number: the Beale number BeN defined in Eq. (1) [14]: 

BeN — f > m/(^swPmean^V) (1) 

The Beale number is almost constant (-0.015) for Stirling en¬ 
gines operating at a temperature ratio z = T c /Th of about 1 /3 and for 
engines which present the same dead volume ratio x and phase 
angle a [14]. 

This relation suggests that the power is related to parameters 
which can be modified to compensate the size effects. As an 
example, the decrease of the power when using smaller displace¬ 
ment (V sw ) can be compensated by an increased frequency or 
pressure. Flowever, their effects for the engine are questionable: 
What are the effects of the frequency and pressure on the 
exchanger or the regenerator performances? Can these compo¬ 
nents be geometrically scaled while maintaining their perfor¬ 
mances? Therefore, the Beale number is incomplete and a coupled 
analysis which includes the engine thermal components is needed 
here. 

Aiming at a more accurate prediction of the power, the Beale 
number has been extended to include the regenerator, the engine 
geometry and the working gas parameters [15,16]. These strategies 
stem on correlations defined from experimental operating mea¬ 
surements of Stirling engines whose displacements range from 
8 dm 3 to 20 dm 3 . The main sources of losses (friction, leakages) can 
then be related to the engine parameters (frequency, pressure, type 
of gas, relative dead volumes). These approaches are accurate as far 
as the derived scaled engines are close to the displacement range. 
The validity of the power prediction and correlation laws for 
significantly smaller volumes may be questionable. Moreover, these 
analyses were not performed for FPSE for which too few experi¬ 
mental data are available in the literature. 

We propose to develop a coupled approach to obtain novel FOM 
allowing easy preliminary designs of FPSE. Each of the parameters 
defined in Fig. 1 will then be evaluated for a given displacement (i.e. 
a given scale) and a given pressure which is usually bounded by 
technological constraints. 

The proposed FOM is based on the similitude theory applied to 
the relevant physics pertaining to FPSE: thermodynamics, conser¬ 
vation of mass, momentum and energy for the gas flow within the 
heat exchangers and the regenerator, mechanical dynamic equi¬ 
librium for the piston and the displacer. Six different types of DGs 
(dimensionless groups) are then obtained. Five DGs are common to 
all types of Stirling engines and one additional DG is proposed to 
account for the dynamic in the case of FPSE. The DGs are evaluated 
for a reference engine. Then, the parameters of the basic compo¬ 
nents of the machine (see Fig. 1 ), can be evaluated keeping constant 
values for the DGs. 


From the derived geometry, the effect of scaling on the classical 
losses in Stirling engine (conduction, regenerator reheat losses and 
pressure drop and gas spring hysteresis) can be studied. The sealing 
of piston is underlined as a major source of losses [2]. In this study, 
this effect is taken into account, evaluating the gas leakage within 
gap between the displacer and the casing (see Fig. 1). Finally, the 
evolution of the power density, the efficiency and the major losses 
are studied. Finally recommendations regarding architecture and 
process are suggested for micro-Stirling engine design. 


2. Similitude approach for the free piston Stirling engines 

Previous works about the similitude theory applied to Stirling 
engines are based on the underlying physics of the Stirling engine 
[17,18]. The basic equations of thermodynamics, gas flow (conserva¬ 
tion of mass, momentum and energy) are written using dimension¬ 
less variables. Dimensionless groups of parameters are then 
enhanced and chosen as representative. To take into account flow 
friction and thermal exchange, usual dimensionless groups are added. 
In addition, the hot and cold temperatures are assumed to remain 
constant during the scaling process and some assumptions on the 
kinematic and the geometry are made. This strategy is recalled here 
and the DGs are established in a comprehensive presentation. 

The obtained DGs were limited to kinematic Stirling engine. In 
this work, the method is extended to the study of the mechanical 
aspects of FPSE and additional DGs are established. 

The proposed dimensionless variables are defined as: 

n * p Ysw 

- Pressure: p = - -= p 

Pmean lvlKl re f 

- Mass within component i : o { = m z /M 

- Mass flow rate through component i : = rhj/(wM) 

-Temperature: T* = T/T ref 

- Time: t* = cot 

- Length: x* = x/L ref 


2.1. Kinematic and geometry 


A kinematic and geometric similarity is assumed. The phase 
angle, the swept volume ratio of the piston and displacer are 
constant. 

Because the ratios of the swept volume to dead volumes are 
supposed to remain fixed, the first DG is then defined as: 




In which i stands for the cooler (i = c), the heater (i = h) or the 
regenerator (i = r) component. Note that only the meaningful terms 
are kept which is why the tt value does not appear in the numerator 
of Eq. (2). 


2.2. Thermodynamics 

From the ideal isothermal model assuming perfect gas law and 
sinusoidal motions, the power output can be expressed. This model, 
used in numerous studies, is called the Schmidt analysis and allows 
analytical expressions for the engine [14]. 

The power can then be expressed as: 


P m — wp m ean^sw 

/ 1 1 VW x(T-l)sin(a:) 

127T2/2 V S w ^/z 2 + k 2 + 2/acos(ce) 



In which t, k and 0 are defined in Eq. (4): 
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T 


K 

P 


T c 

T„ 

^SWC 

^swe 

yj2 tk cos(a:) + t 2 + k 2 
2/l/ S we(Vc + VJc/Tr + V h T) + K + T 



The assumption of kinematic and geometric similarity leads to 
constant k, a and V swe IV sw . Consequently, the Beale number BeN, 
given in Eq. (1), can be defined from Eq. (3). BeN is then chosen as 
the second DG. 


Because the volume can also be expressed as A = V//L/, and 

L ref 3 = V sw , this DG becomes coLi/y/yRT re f 1 /DG V i. The defined DG 
is finally: 



\J T^^ref 


( 11 ) 


Note that DG M . can be seen as a modified Mach number as 
pointed out by Organ [17]. 

The friction factor fy is related to the Reynolds number. From its 
general definition, Re = prUr Hii/M- Using the set of dimensionless 
variables, the Reynolds number can be written as: 


2.3. Gas flow equations 


The energetic similarity is obtained from the ID behavior of an 
ideal gas flow within the heat exchangers and the regenerator. The 
fundamental conservation equations (mass, momentum and en¬ 
ergy) are given by the Eqs. (5)—(7) respectively: 


m+a(V\ = o 

dt dX y-Affi J 



at V%i J a x\A m 7 


dp 

dX 


1 


/f 


771/ 1 771/ 
pAff 2 





+ 



P 

+ —+ 
Pi 



a / 77i/ \ 

ax \A m Ui ) 




In which the free flow area is such as A m = nj 7 rr hi 2 . For the 
momentum Eq. ( 6 ), the fluid dissipation is expressed using the 
Fanning friction factor fy. 

Using the dimensionless variables and noting that for a 
considered ID slug of fluid, one can write: Pi = p/RTi = ( 
\/Aff i )(M/L re f)(da i /dx*) and pfUf = (wM/v 4 ffl )< 7 / Eqs. (5)—(7) are 
written again as: 


-X( d 4] = o 

at* \dx*J ax* 



d&i | d&jUj _ 1 ( A ffiVy RT ref \ L re f 2 dp 


dr ax- 


T 


y toL re f 


f ^ref I'ref ^ l/T 
JF~--— 


Am ax* 2 r hi A ffl 


i \ w i i 


(9) 


a 




dt 




dot 

dX* 




a 

H-tr 

ax* 


v i 

-T + = 

7 2 

/ 




V 




\ 




AmVWFFf 


J 


(T-l) 


ffj 


P 


2 \ 


/ 


\ 


/ 


ffi 




1 _* 7-11 1 

-T + 1 -=5T+~ 

7 7 12 








^ffi \/yRT,ef 


\ 


J 


x (7 - 1) 


T (7/ 

* 

p 


2 \ 




\ 


/ 


" 1-ref > t-r* _ 7* 

p/iCp r hi ! V w 


( 10 ) 


From Eqs. (9) and (10), the DG &jL ref 3 /A ffi \/yRT ve f can be pointed 


out. 


Re = DG M i ai 


2- 1 PmeanUii 


DGyj OJpLf 


( 12 ) 


Because DG Mi and D V i have been already defined, we define: 



PmeanOii 

OJpLf 


(13) 


2.4. Thermal transfer 


The thermal transfer coefficients h in Eq. (10) is related to the 
Stanton number (St = h/pC p u). 

Experimental correlations allow St to be expressed as a 
function of the Reynolds and the Prandtl (Pr = C p p/k) numbers. 
For the heater and the cooler, turbulent flows are reported 
whereas laminar flow usually occurs in the regenerator. Hence, 
the correlations usually used for Stirling engines are chosen 
[13,19] as Eq. (14) for the heat exchangers and Eq. (15) for the 
regenerator. 

St = Pr“ 2 / 3 exp(0.337 - 0.8121og(Re)) (14) 


St = Pr“ 2 / 3 exp( - 0.299 - 0.411og(Re)) (15) 

The ratio h/ppC p L ref /r bi appears in Eq. (10). As a consequence, 
using the previously defined DGs and the Eqs. (14) and (15), two 
new DGs are defined for the heater and cooler (Eq. (16)) and for the 
regenerator (Eq. (17)): 


DG N TUi 



-2/3 



1.812 


^ Pmean ^j 0 812 ,,0.145 




(16) 




DG 


NTUr 



1.412 


fPmean\ 0-412 -0.254 
r 


\ (j) J 


(17) 


2.5. Mechanical dynamics equations 


Dynamic similarity is obtained from the dynamic equilibrium 
equations of the moving parts: Eq. (18) for the piston and Eq. (19) 
for the displacer. 


mpXp = A p p - CpXp + A bp p bp 


(18) 


m d x d 


^rodP Gl^d A b 



f¥m\rhj 

PAt fi 2 


Lf + A vob p bb 


(19) 


Using the dimensionless variables the Eqs. (18) and (19) 
become: 



ApP ref * _ c p • * AbpPref * 
(x) 2 m p L re f wm p p to 2 m p L ref bp 


( 20 ) 
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x* _ ^rodPref * _ , -^d \ /fi ^^ref I I Jj_ r 

d _ w 2 m d L ref P ojm d d m d L ref ^2r hi A ffi 2 p* 1 

^rodPref * 
w 2 m d L ref bd 

( 21 ) 

The natural dynamic of the FPSE is related to the first and the 
last terms of the right hand side of Eqs. (20) and (21). They account 
for the balance between pressure and inertia effects. The effects of 
the dissipative terms are supposed to be negligible for the oper¬ 
ating frequency of the engine. The gas dissipative term has already 
been taken onto account in DG Ri . 

Therefore, the specific dynamic of FPSE can be ensured during 
the scaling process with the following DGs: 


DGd p 

DG Dd 


■^pPref 

9 r 

(22) 

to z m p L re f 

■^rodPref 

(23) 

G)^m d L re f 


3. Method for preliminary design using the DGs 

3.1. Comparison of the DGs for various Stirling engines 

In order to support the proposed approach, the DGs can be eval¬ 
uated for documented engines [17] namely: GPU-3, PD46, V160, 
MP1002CA, USSP40, 400HP. The FPSE RE1000 [20], CTPC [21], BIO 
[22] and DFPSE [23] have also been taken into account. Table 1 pre¬ 
sents some of the relevant characteristics of these engines for scaling. 

The ratio DG Mi /DG vi appears in Eqs. (9) and (10). It is of 
importance in the underlying physics of the engine. Fig. 2a-c 
shows the relation between the DG Mi and DG V i values for the 
regenerator, the cooler and the heater of various real Stirling en¬ 
gines. The different ratios DG Mi /DG V i do not give the same values, 
though a global trend appears from Fig. 2a-c. This analysis backs up 
the proposed approach. 

For the new proposed DGs for FPSE, Fig. 2d presents the ratio 
V^Pmean/rnLref as a function of the frequency which composes the 
dynamic groups DG Dp and DG Dd defined in Eqs. (22) and (23). 
Again, the correlation between these parameters proves that the 
proposed DGs are representative of the FPSE characteristics. 

Similar conclusions can be drawn for the other DGs, but are not 
illustrated here for conciseness reason. 


constant. Relationships between the parameters are obtained from 
the five DGs defined by Eqs. (1), (2), (11), (13), and (17). In logarithmic 
form, the 5 linear equations for 7 unknown variables are: 


-1 -1 -1 

-10 0 

0 0-1 
0 1 -1 


0 -0.412 0.412 


0 0 0 1 

2 110 

0 10 0 

1 -10 0 

-1.412 1.412 0 0 


log(BeN) 

log(DG Vr /4) 

log (dGmt \J yRTt) 

l0g(DG Rer/tt) 

log(DG NTUr /0.74/x -0 ' 254 (Cp/k) - 2/3 ) 


logCVsw) 
log(Pmean) 

log(w) 

l°g( r hr) 

log(Lr) 

log(n r ) 

log(Pm) 


(24) 


The swept volume is representative of the global dimensions of 
the engine and the used pressure, which should be as high as 
possible to increase the performance, can be typically given by 
technological constraints. Therefore, the displacement V sw and the 
pressure p me an are usually fixed. The remaining geometrical vari¬ 
ables U, rhi, Hf, the operating frequency and the estimated power 
can then be evaluated solving Eq. (24). 

For the FPSE, as the operating frequency is given by the previous 
DGs. The values of the moving masses (m p , m d ) and the areas (A p , 
A ro d) have to be evaluated. The logarithmic expressions of the Eqs. 
(22) and (23) are: 


log( m p) = - log(DG D p) + log(p m ean) — 1/3 log(V sw ) 

- 2 log(w) + log (Ap) (25) 


log(m d ) = — log(DG Dd ) + log (p m ean) — 1/3 log(V S w) 

- 2 log(w) + log(A rod ) (26) 

In order to solve these equations, one more assumption is 
needed. The piston and rod surfaces are supposed to scale as U sw 2/3 
(i.e. L re f 2 ). By doing this, the mass values of the piston and the 
displacer are evaluated. This choice may be modified. For example, 
because of technical reason, these surfaces may not be scaled as 
U sw 2/3 leading to other results for the masses. This would be the 
case if pistons are replaced by membranes. 


3.2. Use of the model 3.3. Example of miniature geometries using similitude 

If we consider a given component ( e.g. regenerator) the simili- The scale e is defined from the derivative and the reference 

tude would be achieved during the scaling if the DGs are kept swept volumes as: = V s r ^ f *e 3 . 


Table 1 

Characteristics of the documented Stirling engines. 


Engine 

MP1002CA 

GPU-3 

USSP40 

Philips 400HP 

PD46 

VI60 

RE1000 

CTPC 

B10 

DFPSE 

Type 

Kinematic 






Free piston 




Characteristics DGs 

BeN; DGyj; 

DG Mi ; D^Ri; 

DG NTU j 




BeN; DGyj; 

DG Mi ; DG Ri ; 

DG NTUi ; DG Di 


Gas 

Air 

h 2 

h 2 

He 

He 

He 

He 

He 

He 

He 

Pm (1<W) 

0.25 

8.95 

7.5 

291 

3 

— 

1 

14.5 

— 

— 

Pref (bar) 

14 

70 

150 

110 

102 

150 

70 

150 

1 

10 

/ref (HZ) 

25 

60 

33 

47 

314 

157 

30 

70 

13 

33 

V sw (cm 3 ) 

61.1 

118 

134 

15 300 

77.5 

225 

114.2 

407.6 

22 

384 

L ref (mm) 

39.4 

48.8 

51.1 

248.2 

42.6 

60.8 

48.5 

74.1 

28 

72.7 

m p (kg) if relevant 

— 

— 

— 

— 

— 

— 

6.2 

13.2 

0.5 

6.2 

rrid (kg) if relevant 

— 

— 

— 

— 

— 

— 

0.42 

2.17 

0.08 

0.42 
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DG Mc 



Fig. 2. Energetic correlation for the regenerator of various Stirling engines, a) Regenerator, b) cooler, c) heater, d) dynamic correlation for FPSE. 


For pressures varying from 10 to 70 bars, and displacements 
from 5 to 0.1 mm 3 [e range [0.01, 0.04]), test cases are evaluated 
from the fully documented FPSE RE1000 [20] for the same tem¬ 
peratures (7h ~ 550 °C, T c ~ 60 °C). The scaling strategy consists in 
setting the scale coefficient (i.e. the displacement U sw ) and the 
mean pressure (p me an) and evaluating the operating frequency and 
the geometric parameters for all the heat components (heater, 
cooler, regenerator) and the masses. 

From the detailed engine parameters it is possible to obtain a 
preliminary design of a micro-Stirling defining the lengths, the 


hydraulic radii, the numbers of channels and the masses of the 
engine. Results are given in Table 2. 

The evolutions of the length or hydraulic radius are shown to be 
not proportional to the scale (e.g. comparing cases #1 and #2 in 
Table 2). The regenerator lengths range from 6 to 14 mm for 
pressurized engine, which is consistent with the literature [11,12]. 

As a crucial component, the evolution of the geometry of the 
regenerator is underlined in Fig. 3. Solving Eq. (24) for the DGs values 
of the RE1000, we can estimate that for the regenerator, the hy¬ 
draulic radius rh r scales as Pmean -0 58 ^sw 014 (orp me an _ 0 ' 58 c 0 42 


Table 2 

Geometric and operating parameters for four derivated micro-engines. The bold figures are the chosen parameters from which the similitude process is performed. 



Reference (RE1000 [20]) 

Scaling case #1 

Scaling case #2 

Scaling case #3 

Scaling case #4 

Swept volume V syv [mm 3 ] 

69 000 

5 

0.5 

0.5 

0.1 

Scale e 

1 

1/24 (— 0.04) 

1/52 (— 0.02) 

1/52 (— 0.02) 

1/88 (— 0.01) 

Op. pressure p mean [bar] 

70 

70 

70 

10 

10 

Op. frequency/= w/2tu [Hz] 

30 

191 

300 

134 

183 

Displacer 0 Dd [mm] 

57.2 

2.3 

1.1 

1.1 

0.6 

Cooler 

Length L c [mm] 

79.2 

10.2 

5.24 

12.4 

8.8 

Hyd. radius r^ c [pm] 

223 

72 

55 

193 

160 

Heater 

Length Lh [mm] 

183.4 

23.7 

16.7 

20.4 

14.4 

Hyd. radius rhh [pm] 

590 

191 

158 

720 

595 

Regenerator 

Length L r [mm] 

64.4 

10 

6.4 

14.4 

10.5 

Hyd. radius [pm] 

70 

19 

14 

43 

34 

Wire 0 d w [pm] 

89 

24 

17 

54 

43 

Porosity If [%] 

75.9 

75.9 

75.9 

75.9 

75.9 

Max. Knudsen number (regenerator) 

8.63 10~ 5 

3.2 10~ 4 

4.4 10~ 4 

9.9 10~ 4 

1.23 10 3 
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using the scale e) and the length L r as Pmean“ 0 ' 41 ^sw° 19 or 
Pmean e ° 57 * The test cases #1-4 are presented as black circles. 

The black line is the usual proportional scaling law for 70 bars 
(the length scales as e). According to the similitude approach, the 
proposed dimensions theoretically ensure a similar effectiveness 
compared to the reference engine. Fig. 3 clearly shows that pro¬ 
portional scaling would not lead to efficient regenerator. Moreover, 
moderate pressures tend to increase the dimensions. 

3.4. Model assumption validity 

The proposed approach relies on the consistency of the 
continuous flow theory. The minimal dimensions of the engine that 
can be designed using the proposed FOM are then defined here as 
the lowest swept volume associated to the rarefaction effect limit. 

It can be assessed evaluating the Knudsen number Kn. The 
minimal equivalent hydraulic radius is the reference length for the 
regenerator component. As an example, it has been given for the 
four test cases in Table 2. If the mean free path of the molecule is far 
smaller than the characteristic length of the considered problem 
(Kn < 1CT 3 ), then the continuum theory can be applied and so the 
proposed model. The Knudsen numbers are evaluated for the 
regenerator and are given in Fig. 4. 

It appears that a size reduction of about e = 0.01 (U sw = 0.1 mm 3 ) 
for a minimal operating pressure of 20 bars, represents the validity 
limit for the model. Under these dimensions, usually neglected 
phenomena may have to be taken into account [24]. 

Some theoretical works studied the effect of quantum de¬ 
generacy of the performances of Stirling cooler. Quantum de¬ 
generacy of gas appears at very low temperatures or for dimensions 
under micrometer scale. Although, some results suggest that 
operating a Stirling cycle machines at nanoscale is possible [25,26], 
such studies are beyond the scope of this paper. 

4. Scaling laws for FPSE 

Scaling analysis for Stirling machines has been proposed in Refs. 
[10-12]. They are based on the effects of a geometrical scaling on 
the various relevant phenomena for the Stirling machine as 
described in Table 3. In these studies, any length scale as e, area as e 2 
and volume as e 3 . The scaling laws derived from the expression of 
the different effects and losses are called “CSL” (classical scaling 
laws) in the following. 

Using the proposed FOM which takes into account the coupling 
in Stirling engine, NSL (new scaling laws) are defined. They stem on 



0.01 1 100 10 4 

V sw (mm 3 ) 

Fig. 4. Evaluation of the Knudsen number associated to the regenerator. 

the introduction of the solution of Eqs. (24)—(26) in the expression 
of the losses and power. 

In this study, the effect of scaling on the frequency, the regen¬ 
erator losses and the output power are taken into account. Addi¬ 
tionally, the effect of a gap between the displacer and the housing is 
considered. This is thought to be especially relevant because of the 
problem of sealing and manufacturing tolerance at small scale. The 
output power and losses are also proposed as specific quantities by 
dividing their expression by the representative volume V sw . 

4.1. Operating frequency and specific power 

The operating frequency is related to the square root of the 
stiffness to the mass ratio. Whether associated to the elastic strain 
of material (mechanical spring) or the adiabatic gas-spring, the 
stiffness scales as e [10]. Because, for CSL, the mass scales as 3 [10], 
the frequency law is e~\ This evolution is the black curve plotted in 
Fig. 5a. 

Using the FOM, the natural frequency w is given by the solution 
of the set of Eq. (24) whenp mea n and e (e = l/ sw 1/3 ) are fixed, w is a 
then function of the DGs, the gas characteristics, p me an and e. Then, 
the new scaling law for the frequency is calculated as: 

(V — P mean £ 

The evolutions of the frequency according to the NSL are plotted 
in Fig. 5 for different pressures. Compared to the CSL, the operating 
frequency which theoretically ensures the same level of 



0.1 1 10 100 1000 10 4 

V s „ (mm 3 ) 

Fig. 3. a) Length of the regenerator as a function of the displacement Vsw, 



0.1 1 10 100 1000 10 4 

Vsw (mm 3 ) 

hydraulic radius for the regenerator as a function of the displacement V sw . 
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Table 3 

Scale effects on the various relevant physical phenomena for Stirling machines. 
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Phenomena 


Dynamic parameters 
Natural frequency 

Mass 

Regenerator 
Gas flow 

Pressure drop 

Regenerator reheat thermal loss 

Conduction 

Gap 

Gas flow 

Leakage reheat thermal loss 
Gas-spring 
Hysteresis gas loss 
Power 

Output power 
Specific quantities 
Specific output power 
Hysteresis gas loss 

Conduction loss 

Regenerator reheat thermal loss 
Leakage reheat thermal loss 
Viscous dissipations 


Expression Units Classical Scaling law Comments 

scaling law from FOM 


[K 

6J V M 

rad s 1 

[10] 

n 0.41 ,,-0.58 

Pmean 


in p = pV p 

kg 

e 3 [10] 

n 0.18 „2.16 
Pmean £ 

For piston or displacer 

^^PmeanKef [n| 

Ki ref 

m 3 s" 1 

Pmean^ 

n 1.41 2.42 

Pmean £ 

Constant T re f 

W^[12] 

pAf r h 

Pa 

Pmean ^ 

n 2.16 „0.18 
Pmean £ 

Constant/f 

Qr = P -tj r )mC p AT[12] 

W 

PmeanC 

n 1.41 2.42 

Pmean £ 

Constant AT % is the 

FA,- 




regenerator efficiency 

Qcond =-j—kT [10—12] 

W 

e 

n 0.83 1.83 

Pmean £ 

Cst AT For the regenerator 

L r 




characteristics 


3 -1 

m s 

Pmean 

n 2.57 „1.59 

Pmean £ 

Gap leakage assuming that 
the gap scales as e 1/3 

Qgap = rnCpAT 

W 

Pmean^£ 

n 2.57 „1.59 

Pmean £ 

For gap leakage, with cst AT 

Qhys — \J 22 ^7 (^ UT/vPmean^^ ^ [14] 

W 

n 15 

V mean^ 

n 0.71 1.71 

Pmean £ 


Pm 

W 

Pmean^ 

n 1.41 2.41 

Pmean £ 

From BeN 

Pm = Pm/Vsw 

W m 3 

PmeanC 

n 1.41 —0.59 

Pmean £ 


Qhys = Qhys/^sw 

W m 3 

PmeanC 

n 0.71 —1.29 

Pmean £ 


Qcond = Qcond/^sw 

W m 3 

-2 

£ 

n 0.83—1.17 

Pmean £ 

Constant AT 

Qr* = (l -p r )mC p AT/V sw 

W m 3 

PmeanC 

n 1.41 —0.58 

Pmean £ 

With constant AT and ?] r 

Q^mCpAT/Vsw 

W m 3 

Pmean^ 

n 2.57-1.41 

Pmean £ 

For gap leakage, with cst AT 

Kb=J*P/V sw 

W m 3 

Pmean 3 £ ^ 

n 3.57—0.4 

Pmean £ 

For the regenerator flow, 
with scaled Ap 


W m 3 

Pmean^e ^ 

n 4.73-1.23 

Pmean £ 

For the gap leakage flow, 
with scaled Ap 


performances for the exchangers (i.e. heat transfer coefficient and 
friction factor) has to be much lower. Therefore, the mass of the 
pistons will have to be chosen properly which is the point of DG Dp 
and DGod and will be discussed in the following. 

The operating frequency of the FPSE partly relies on the moving 
masses. Fig. 6 shows the evolution of the displacer mass with 
respect to the scale e. The thick black line stands for the CSL (e 3 as 
can be read in Table 3). Using the dynamical results from similitude, 
the mass scales as Pmean 018 c 2 14 . For small scale variation, the 
evolutions are close (Fig. 6a). However, for the smallest engines, 
drastic differences can be seen in Fig. 6b. Assuming that the stiff¬ 
ness effect is associated to gas springs, the proposed new scaling of 
the displacer mass ensures the operating frequency to be optimal. It 



Fig. 5. Operating frequency for pressure from 10 to 70 bars. 


can also be inferred that the relative mass with respect to the en¬ 
gine size will be more important in a micro-Stirling than in its 
macro-scale counterpart. 

The output power is related to the operating frequency. Fig. 7 
presents the CSL (black curve) and the NSL for the specific output 
power. The power density increase when miniaturizing is predicted 
by both the scaling laws. However, on the one hand, the frequency 
has to be tuned; on the other hand the geometrical parameters 
have to vary properly to maintain the heat exchange properties. 
This is ensured by the proposed laws whereas the classical law 
would lead to too high frequency (Fig. 5) and too small geometry 
(Fig. 3). 

Therefore, the proposed evolution of the power density is 
thought to be more relevant. Still, cubic millimeter displacement 
engine (e ~ 0.04) can present up to 7 times more specific power 
than macro-scale engine. The role of the pressure is presented in 
Fig. 7. The power density is shown to be not proportional to the 
pressure which is a new result. Again, the maximal attainable 
pressure is to be sought to increase the specific power of the 
engine. 

4.2. Losses 

The FOM theoretically ensures the same behavior for the heat 
exchangers and the regenerator. Then, the effectiveness of the 
regenerator, the heat exchange coefficient and the friction factor 
are theoretically constant. However, the losses phenomena have to 
be analyzed to anticipate the engine efficiency. Using a previously 
developed model for FPSE [9], the losses have been evaluated for 
the RE1000 and are used here to obtain their initial distributions at 
full scale [e = 1, subscript FS in the following). Then, their evolu¬ 
tions with scaling can be anticipated using the NSL. 
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030 0.50 0.70 1.00 1.50 0.015 0.020 0.030 

8 £ 

Fig. 6. a) Displacer mass as a function of the scale for e range [1/4, 3/2], b) for e range [1/88,1/24]. 



Fig. 7. Power density evolution for pressure from 10 to 70 bars. 


4.2.1. Conduction losses 

The CSL for the heat conduction across the regenerator suggests 
that the associated specific heat loss scales as e~ 2 as can be read in 
Table 3 (black curve in Fig. 8a). Following the similitude, it can be 
inferred that the same loss scales asp mea n° 83 c _117 .The conduction 
loss increase is smaller than for the CSL because the regenerator 
length scales as e 0,57 , (see Figs. 3 and 8a). Additionally, this new 


result underlines the indirect effect of the pressure which could be 
used to reduce the conduction loss. 

Fig. 8b plots the evolution of the percentage of the specific 
conduction losses relative to the total losses. It can be seen that the 
conduction can account for up to 60% of the total losses for mm 3 
displacement engine (e ~ 0.04) and is superseded by other losses 
for smaller engines. Note that the conduction loss influence is 
enhanced by low pressure which means that other losses may 
drastically lessen with pressure. 

4.2.2. Regenerator losses 

Regenerator losses are: the reheat thermal losses associated to 
its effectiveness and the pressure drop. Its effectiveness % is 
theoretically constant since it can be expressed as a function of the 
NTU [12] which is taken into account in DGntup 

The increase of the thermal loss when miniaturizing (Fig. 9a) 
appears to be moderate in comparison with the evolution of the 
conduction loss. Moreover, from the full scale value, the weight of 
this loss is only 0.7% of the total loss and this value decreases as 
other losses prevail when miniaturizing (Fig. 9b). 

The NSL for the pressure drop loss is obtained assuming a 
constant friction factor coefficient which is in accordance with the 
constant DG Rr . The increase of this loss appears to be much less 
than for the CSL (Fig. 9c). Moreover, reducing the pressure can 
drastically diminish this loss. Using a pressure of 10 bars makes the 
contribution of the pressure drop loss in the regenerator less than 
1% of the total losses (Fig. 9d). 
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Fig. 8. a) Ratio of the specific conduction loss evolution to the full scale specific conduction, b) percentage of the specific conduction losses relative to the total losses. 
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Fig. 9. a) Relative specific regenerator reheat thermal loss evolution, b) evolution of the ratio of the reheat thermal loss to the total losses, c) relative specific pressure drop loss 
evolution, d) evolution of the ratio of the specific pressure drop loss to the total losses. 


4.2.3. Gas spring hysteresis 

Another damping effect is the gas spring hysteresis. For FPSE 
it can also modify the phase angle and the amplitudes of the 
piston and displacer, which affects both power and efficiency. 
From Fig. 10a, it can be seen that the miniaturization greatly 
increases the hysteresis loss. Flowever, in the case of the RE1000, 
this loss is small compared to the total of losses as can be seen 
in Fig. 10b for e = 1. As long as the RE1000 is used as a 


reference engine for preliminary design, the gas-spring loss is 
not critical. 

4.2.4. Effect of the gap gas leakage 

From a technological point of view, materials and fabrication 
constraints are challenging when designing micro-engines. From the 
theoretical prospective studies, friction losses and leakage through 
the cylinder-piston gap was underlined as one of the main 
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Fig. 10. a) Specific gas spring hysteresis losses evolution, b) evolution of the ratio of the specific hysteresis losses to the total losses. 
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0.02 0.05 0.10 0.20 0.50 1.00 


e 

Fig. 11. a) Relative specific gap reheat thermal loss evolution, b) evolution of the ratio of the 
the gap, d) evolution of the ratio of the specific pressure drop loss to the total losses. 


limitations [1,2,27,28]. Gas leakage through the displacer-casing 
clearance induces reheat losses due to the flow of cold gas into the 
hot expansion chamber which impact the efficiency of the engine. In 
classical scale analyses, the geometrical gap is supposed to scale as e. 
We choose here to highlight the technological point of view. We 
assume that no sealing component is used. From the engineering 
tolerance theory [29], a classical IT5 quality fitting is supposed to be 
adopted and it can be inferred that it scales as e 1/3 . 




0.02 0.05 0.10 0.20 0.50 1.00 
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8 

reheat thermal loss to the total losses, c) relative specific pressure drop loss evolution for 


Due to the gap between the displacer and the casing a gas flow 
occurs. Its evaluation can be obtained using the hydrodynamic 
lubrication theory [30]. In addition to the pressure drop, a reheat 
thermal loss is considered. The gas flow through the gap strongly 
increases when miniaturizing and both the reheat and pressure 
losses increase (Fig. 11a and c). 

The reheat thermal loss has the most important detrimental 
effect as can be seen in Fig. lib. If the pressure is decreased, 
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Fig. 12. a) Relative efficiency, b) relative specific power as a function of the relative efficiency for a scale ee [0.01,1]. 
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the effect of the gap weakens. In this case, the conduction loss 
becomes the major source of loss (see Fig. 8b forp me an = 10 bars). 


4.3. Efficiency 


The efficiency of the micro-Stirling engine is estimated here as 
the output mechanical power P m divided by the total input thermal 
power following Eq. (27): 


V 


Tn 


(Pm + P\ 


(27) 




in which 7)c is the Carnot efficiency. 

The relative efficiency with respect to the scale e is shown in 
Fig. 12a. Using CSL would lead to small efficiency (black curve in 
Fig. 12a). The study of the various losses showed that a moderate 
pressure tends to reduce the increase for all the losses and a better 
efficiency can be obtained. 

Fig. 12b presents the specific power and efficiency evolution 
with respect to e. Though, the efficiency is improved with lower 
pressure, the effect is much more important on the specific power 
and reduces it drastically. As an example forp mea n = 40 bars (doted 
red curve in the web version), the efficiency is 30% lower than for 
Pmean = 10 bars but the power is six times greater. 

As a conclusion, some hints for designing micro-Stirling engine can 
be deduced form this analysis. First the leakage between the displacer 
and the housing induces the main source of loss. Sliding piston should 
be avoided. Low pressures reduce this loss but have a great impact on 
the power. The conduction loss appears to be the second most 
important source of loss. Therefore, if membranes are used instead of 
piston, thermal conduction through the membrane thickness would 
result in much higher conduction losses. It appears that pressurization 
is more important in micro-Stirling than in its macro-scale counter¬ 
part. It may be an explanation for the fails of the few experimental 
micro-Stirling design attempts using atmospheric pressure [31,32]. 

From the dynamical point of view, it was demonstrated that the 
inertia effect has to be stronger for miniature engines than for their 
macro-scale counterparts. New membrane structures are then to be 
designed to match the thermal isolation and the inertia re¬ 
quirements, such as proposed in Ref. [33]. 


5. Conclusions 


The effects of size reduction on performances have been studied 
for the FPSE, which is the best suited Stirling architecture for 
miniaturization. From the similitude approach a simple first order 
model has been proposed. It allows a preliminary design to be 
obtained from a full scale engine. The geometrical parameters 
(length, hydraulic radius, number of channel) for the heat ex¬ 
changers and the regenerator as well as the operating frequency 
and the masses of the pistons are obtained. 

Moreover, the evolution laws of these parameters allow new 
scaling laws for the power and the major sources of losses to be 
defined. Conclusions from these laws suggest that the specific power 
of a Stirling engine can be significantly increased with size reduction 
and using high pressure. However, the increase is less than predicted 
by the classical geometrical scaling law and the new results are 
thought to be more realistic. Besides, the role of the pressure is 
underlined and clearly stated by the proposed scaling laws. Matching 
the optimal operating frequency allows to maintain high effective¬ 
ness for the heat exchangers and the regenerator. It requires that the 
moving masses are properly designed. The relative inertia with 
respect to the engine dimensions increases as the size decreases 
which lead to larger moving components. 


The detrimental effect of gap leakage on the engine efficiency 
has been underlined. Because the gap induces most of the losses, 
the sliding piston architecture appears to be not suitable for micro¬ 
engine. The use of membrane instead of the piston and the dis¬ 
placer could be a solution to this issue. However, the membrane 
would have to ensure high thermal isolation and to match the 
frequency requirement which means high inertia effect. 
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